m 



PROCEEDINGS 

°F SCIENCE 



PoS (Confinement X) 111 



Determination of SU(2) ChPT LECs 

from 2+1 flavor staggered lattice simulations 



Enno E. Scholz", Szabolcs Borsanyi^, Stephan Diirr'^' , Zoltan Fodor* '^'', 
Stefan Krieg^ Andreas Schafer ', and Kalman K. Szabo^ 

"Universitat Regensburg, Universitdtsstr. 31, D-93053 Regensburg, Germany 
^Bergische Universitat Wuppertal, Gaufistr 20, D-42119 Wuppertal, Germany 
'^Julich Supercomputing Centre, Forschungszentrum Jillich, D-52425 Jiilich, Germany 
Institute for Theoretical Physics, Eotvos University, H-1 117 Budapest, Hungary 



E-mail: 


snno . scholz (AT 


physik . uni-regensburg . de 


durr (AT) itp . unibe . ch 





By fitting pion masses and decay constants from 2+1 flavor staggered lattice simulations to the 
predictions of NLO and NNLO SU(2) chiral perturbation theory we determine the low-energy 
constants and £4. The lattice ensembles were generated by the Wuppertal-Budapest collab- 
oration and cover pion masses in the range of 135 to 435 MeV and lattice scales between 0.7 
and 2.0 GeV. By choosing a suitable scaling trajectory, we were able to demonstrate that precise 
and stable results for the LECs can be obtained from continuum ChPT to NLO. The pion masses 
available in this work also allow us to study the applicabiUty of using ChPT to extrapolate from 
higher masses to the physical pion mass. 



Xth Quark Confinement and the Hadron Spectrum, 
October 8-12, 2012 

TUM Campus Garching, Munich, Germany 
* Speaker. 



© Copyright owned by the author(.s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Determination ofSU(2) ChPT LECs from 2+1 flavor staggered lattice simulations 



Enno E. Scholz 



1. Introduction 

Chiial perturbation theory (ChPT) [|l} ^ is a widely used tool in many phenomenological 
applications and also helpful to guide an extrapolation to lighter quark masses in lattice-QCD 
simulations. Here we will report on a determination of the NLO low-energy constants (LECs) £3 
and £4 which appear in the light quark mass dependence of the pseudo-scalar meson masses and 
decay constants in SU(2) ChPT. 

We analyze configurations generated by the Wuppertal-Budapest Collaboration [Q, ^ ^ ^ ^ ^ 
using the Symanzik glue and 2-fold stout-smeared staggered fermion action for a 2-i- 1 flavor QCD- 
simulation. The mass of the single flavor has been kept at the value of the physical strange quark 
mass, whereas the two degenerate lighter quark masses have been varied such that light meson 
masses in the range of 135 to 435 MeV were simulated. The simulations were performed at six 
different gauge couplings jS, resulting in lattice scales between 0.7 and 2.0 GeV (see next section 
for details on how the scale has been determined). Figure |l] shows a landscape plot of our simulated 
pion masses squared versus the lattice spacing. Details about the simulated gauge couplings, lattice 
volumes, and tuning of the input quark mass values are reported in our publication |^. 

The 2-fold stout-smeared version of the staggered quark action has been proven to be advan- 
tageous [^] in reducing the inevitable taste-breaking of staggered femiion formulations. Therefore, 
in this work we only consider the pseudo-scalar mesons with taste matrix 75 when measuring me- 
son masses or decay constants. Again, for details of the computation of these quantities we refer 
the reader to [Q. 

2. Scale setting and physical quark masses 

To set the scale at each simulated gauge coupling j8 and identify the physical point, i.e. the 
average up/down quark mass m^^^^ = (mu + md)/2 corresponding to a pion in the isospin limit 
with an estimated mass of M^^^ = 134.8 MeV [[Io|], we use a two-step procedure. First, we ex- 
trapolate the ratio {aMji)'^ / [afj^Y of the measured squared meson masses and decay constants 
to its physical value (134.8MeV/130.41 MeV)^ = 1.06846, where we also used the PDG-value 
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Figure 1: Simulated pion masses squared 
vs. lattice spacing a at six different 
gauge couplings j3. Horizontal dashed 
lines indicate cuts on the mass range in our 
ChPT fits. The lower left corner marked 
by blue solid lines is our final preferred fit 
range. See text for details. 
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Figure 2: Left panel: ratio {aMj^)^ / {af„)^ extrapolated to {Alt^' /f^""^')^ = 1 .06846 to obtain amP'^y', right 
panel: afj^ extrapolated to am^^^^ to obtain l/a; both at j3 = 3.85. Points marked by black symbols are 
included in the fit, while those marked by green symbols are excluded. 



yphys _ 230.41 MeV [11 1. In that way am^^^^ is obtained. In the second step, we extrapolate afji 
to this quark mass value and obtain the lattice scale with the help of the PDG- value f^"^^ . For the 
extrapolation we used two different ansatze: a quadratic and a rational (linear in numerator and 
denominator) fit form. An example of these extrapolations is shown for the ensembles at j8 = 3.85 
in Fig. |2[ There, only the five lightest points were used in the fits, which is a typical choice for 
all other ensembles. We stress that here, like in the chiral fits to be discussed below, the data has 
been corrected for finite volume effects beforehand, by means of using the two- and three-loop re- 



summed formulae of [ |12| ] for the pion decay constants and masses, respectively. Our spatial lattice 
volumes L are in the range (4.3fm)3 - (6.8fm)3 with a minimal M^L ^ 3.3. This ensures that we 
only observe small finite volume corrections. In case of the pion mass the correction factors vary 
between 0.1 and 2.7 per-mille and in case of the decay constant between 0.2 and 7.5 per-mille. 

By fixing l/a and am^^^^ in the way described above, the meson masses and decay constants 
show no discretization effects at all directly at the physical point and we can assume those effects 
to be small (since of higher order in the quark masses and/or lattice spacing) in the vicinity of the 
physical point, i.e. in the mass range covered by our fits. Such discretization effects, of course, are 
present in other obsei-vables, which are not considered in this work. 



3. Fits to NLO SU(2) ChPT 

The quark mass dependence of the finite-volume corrected data for the meson masses and 
decay constants is fitted simultaneously at different jS -values using the NL0-SU(2) ChPT formulae 
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where we made use of the already determined Xja and ayrP'^'^^ to scale the quark masses and 
the meson masses and decay constants measured in lattice units. This fit has four free param- 
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Figure 3: Combined global unconstrained NLO fit of the decay constants (left panel) and squared meson 
masses divided by the quark mass ratio (right panel) to ensembles with l/a> 1.6 GeV and 1 35 MeV < M,i < 
240 Me V. Points included in the fit range are marked by black symbols, points excluded by green symbols. 



eters: two NLO low-energy scales A3, A4 (related to the LECs ii = log[A?/(Mr^')2]), the de- 
cay constant in the SU(2) chiral limit / and the renormalization scheme-independent combination 
^phys _ (2BmP^'^^) of the LO low-energy constant B and the physical quark mass m^^'^^. Since we 
used the physical pion mass and decay constant to set the scale at each set of ensembles with a 
common j8 , and furthermore each set contains at least one data point in close vicinity to the phys- 
ical point, our ansatz should reproduce the physical point. Therefore, we also used a parameter- 
reduced chiral fit, where this constraint has been implemented and only two free fit parameters 
remain, which we chose to be ^P^y^ and /. Our exact implementation of the parameter-reduced fit 
formulae is reported in [^. 

We would like to point out that the chiral fit formulae do not include any taste breaking effects, 
i.e., we did not use staggered ChPT. This seems justified to us, since we are only considering 75- 
taste mesons as mentioned above and use these to define our scaling trajectory at the physical 
point. In other words, since the meson mass and decay constant at the physical point were used to 
set the quark masses and lattice scales, no discretization or taste breaking effects are present in the 
ChPT formulae for and fji as discussed above. Furthermore, taste breaking effects are reduced 
anyway by the choice of the fermion action as mentioned above. 

For the chiral fits, we first applied several different cuts for the heaviest mass included in the fit 
range. These cuts are indicated by the horizontal dashed lines in the landscape plot of our simulated 
points. Fig. [l[ In a next step, we studied the effects of excluding one or more lattice spacings from 
the fitted data. In the end, judging by the resulting (uncorrelated) ;^^/d.o.f of the fit and reaching 
a plateau for the fitted parameters, we chose 135MeV < Mj^ < 240 MeV and l/a > 1.6GeV to 
be our preferred fit range. This choice is indicated by the solid blue lines in Fig. |l|. The resulting 
combined global fit to the unconstrained fit formulae is shown in Fig. |3[ In Figure ^ we show the 
impact of the different mass cuts on the fitted parameters and the resulting (uncoiTelated) ;^^/d.o.f. 
using only ensembles with l/a > 1.6 GeV. For the systematic uncertainty we decided to take the 



4 



Determination ofSU(2) ChPT LECs from 2+1 flavor staggered lattice simulations 



Enno E. Scholz 



range [MeV] - ^ 
[ 135 : 160] - 
[ 135 : 195] - 
[135:240] - 
[135:275] - 
[135:340] - 
[135:390] - 

[135:435] - 

0.0182 



■ [GeV^ 



1/a> 1.6GeV 



range [MeV] 
[ 135 : 160] 
[ 135 : 195] 
[135 :240] 
[135 :275] 
[135 :340] 
[135 :390] 
[ 135:435] 



f[IVIeVl 1/a 


>1.6GeV 




























IK 




* 







0.0184 0.0186 O.Oie 



120 



range [MeV] 


[135 


160] 


[ 135 


195] 


[135 


240] 


[135 


275] 


[ 135 


340] 


[135 


390] 


[135 


435] 



A3 [l^eV] 



1/a > 1.6 GeV 
^* 



400 450 500 550 600 650 700 750 800 



M^ range [MeV] - 


[135 


160] - 


[135 


195] - 


[135 


240] - 


[135 


275] - 


[135 


340] - 


[135 


390] - 


[135 


435] - 



1/a > 1.6 GeV 



900 



1000 



1100 



1200 



M, range [MeV] - x 'd.o.f. 

[135 : 160] ^ 

[ 135 : 195] — JK- 

[ 135:240] - -^K— 

[135:275] - ■ X ■ 

[135:340] - ' X ' 

[ 135 :390] - -^K- 

[135 :'*35] [ 



1/a > 1.6 GeV 



0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 



Figure 4: Dependence of the fitted parameters and (uncor- 
related) j^/d.o.f. on the mass range of fits only including 
ensembles with 1/a > 1.6 GeV. The solid, dashed, and 
dashed-dotted blue lines show our central value, statistical 
and statistical plus systematic error bands, respectively, as 
determined from these fits. See text for details. 



variation of all these points with respect to the central value from our preferred fit. The central 
value and statistical uncertainty from our preferred fit is indicated in each panel by a solid and 
dashed line, respectively, whereas the total (statistical plus systematic) uncertainty is indicated by 
a dashed-dotted line. Details on the impact of removing one or more lattice spacing from the fits 
can be found in [^. 

We repeated the same analysis with the parameter-reduced SU(2) ChPT-formulae, which have 
been constrained to hit the physical point. In Table [l| we summarize the obtained fit parameters from 
each ansatz. Note that for the parameter-reduced fit only the first two lines are fitted parameters 
while the two low-energy scales A3, A4 are derived from these. The table also contains our final 
estimate for the NLO SU(2) LECs, which were obtained by combining the results from the two fit 
ansatze in the following way: we averaged the central values and the statistical uncertainties. For 
the square of the systematic uncertainty we sum the squares of the average systematic uncertainty 
and the spread of the central values. 



unconstrained parameter-reduced 



final 



;(.phy.s/(io-2GeV 
//MeV 
A3 /MeV 
A4/MeV 



1.8578(17)(39) 
122.70(08) (41) 
628(23) (57) 
1,012(16)(83) 



1.8639(18)(44) 
122.73(06) (28) 
678(40)(119) 
1,006(15)(71) 



1.8609(18)(74) 
122.72(07) (35) 
653(32)(101) 
1,009(16)(77) 



Table 1: Results with statistical and systematic uncertainties for fitted parameters from unconstrained (first 
column) and parameter-reduced (second column) NLO SU(2) ChPT fits. Note that in the parameter-reduced 
case only the first two parameters are actual free fit parameters, while the other two are derived therefrom. 
The third column shows our combined final results, see text for details. 
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Figure 5: Results for the LECs (left panel) and £4 (right panel) from unconstrained NLO SU(2) ChPT 
fits for different mass ranges including (above top-most dashed horizontal line) and excluding near-physical 
masses. Only ensembles with l/a> 1.6GeV were included in the fits. Blue lines indicate our central value 
and error bands in this set-up. 
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Figure 6: Extrapolated physical pion mass squared (left panel) and decay constant (right panel) from uncon- 
strained NLO SU(2) ChPT for different mass ranges including (above top-most dashed horizontal line) and 
excluding near-physical masses. Blue lines indicate the experimental values and error bands from [|ol |ll]|. 



In addition to our main NLO SU(2) ChPT fits, we also varied the lower cut on the pion masses 
in the unconstrained fits. Especially, this is of interest because nowadays many lattice simulations 
still do not include the physical point and use ChPT to extrapolate towards it. In Figure || we 
show how, e.g., the low-energy constants £3 and £4 (related to the scales A3 and A4, respectively) 
change due to the various lower and higher mass cuts (in the upper part of each plot (grey shaded 
area) the results of fits including the physical point are shown, which have been used to obtain the 
systematic uncertainty, see above). One can see that £3, which predominantly influences the quark 
mass dependence of M^, is still in agreement with results from fits including the near physical 
points, whereas £4, which predominantly influences the quark mass dependence of f^, shows some 
deviations when the physical point is excluded. In Figure ^ we compare the pion mass and decay 
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Figure 7: Combined global fit to NNLO SU(2) ChPT using priors for A12, A:^2 and kf. Only ensembles 
with l/a> 1.6 GeV and 135 MeV < Mj^ < 340 Me V are included in the fit (only data points included in the 
fit range are shown in these plots). 
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Figure 8: Results for A3 (left panel) and A4 (right panel) from NNLO SU(2) ChPT fits to different mass 
ranges using priors for A12, kj^i, and kf. Only ensembles with l/a > 1.6 GeV were included in the fits. Blue 
lines indicate our final central value and error bands. 



constant extrapolated from our unconstrained NLO SU(2) ChPT fits (including and excluding the 
near physical points) to the experimental values [|^, 11]. Also here the pion decay constant shows 
more deviations from the expected result, once more and more lighter masses are excluded from 
the fits. In our opinion, these observations illustrate the danger inherent in applying NLO SU(2) 
ChPT-formulae to lattice data lacking data points with light enough pion masses. 



4. Fits to NNLO SU(2) ChPT 

To check the influence of higher orders in SU(2) ChPT on our results for the LECs, we also 
extended our fits to include the next-to-next-to-leading order (NNLO) contributions. In our set-up 
for the fits of the meson masses and decay constants, we had to add three new fit parameters: the 
low-energy scale A12, which we defined as a combination of the NLO low-energy scales Ai and A2: 
logAj2 = (71ogAj + 8 log A2)/ 15 and two parameters for the NNLO-contributions in the meson 
mass and decay constant dependence, kj^i and kf, respectively (see for details). Unfortunately, 
our amount of data at fine enough lattice spacings and light quark masses did not turn out to lead 
to stable fits. Therefore, in the end we opted for using priors on the three additional fit parameters 
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in the NNLO-fits. Details on the choice for the priors are reported in In Figure ^ we show a 
typical fit of our data to the NNLO formulae using priors. Since we were not able to resolve the 
new parameters in a satisfactory way without a priori input we refrain from quoting results for these 
pai^ameters. But we would like to stress that the change in the determined values for the NLO-LECs 
is only minor as can be seen from Fig. || where we show the results for A3 and A4 from a NNLO-fit 
(using priors for kj^^i, kf, and A12) at different mass ranges. There the vertical Unes show our final 
results and error bands from our NLO-fits. 

5. Conclusions 

We determined the following SU(2) LECs from our NLO ChPT fits to meson masses and 
decay constants measured on staggered 2+1 flavor lattice simulations of QCD, cf. also the right- 
most column of Tab. |l|. 

IBm^^y' = 1.8609(18)s,at(74)syst • 10-2GeV2, / = 122.72(07)stat(35)systMeV, 

£3 = 3.16(10)stat(29).syst, ^"4 = 4.03(03)stat(16)syst. 

(Here we quoted the LECs £t, and £4, which are related to the scales A3 and A4, respectively.) In 
addition, we also provide the ratio of the extrapolated physical decay constant to its value in the 
chiral limit f^^^ / f = 1.0627 (06 )stat (27 )syst- Using the value determined for the average light 
quark mass in [ |T3| , pl| ], we obtain from our fitted value for j(^^y^ the condensate parameter 

gMS,2GeV ^ 2.682(36)stat(39)systGeV, I^.2GeV ^ (272.3(1.2),tat(1.4)systMeV)', 

where the condensate was obtained by multiplying B with our result for /^/2. 

The speaker acknowledges support from the DFG SFB/TR 55 and the EU grants PIRG07-GA- 
2010-268367 and PITN-GA-2009-238353 (ITN STRONGnet). 
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